Unusual Response to a Localized Perturbation in a Generalized Elastic Model 
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The generalized elastic model encompasses several physical systems such as polymers, membranes, 
single file systems, fluctuating surfaces and rough interfaces. We consider the case of an applied 
localized potential, namely an external force acting only on a single (tagged) probe, leaving the rest 
of the system unaffected. We derive the fractional Langevin equation for the tagged probe, as well 
as for a generic (untagged) probe, where the force is not directly applied. Within the framework of 
the fluctuation-dissipation relations, we discuss the unexpected physical scenarios arising when the 
force is constant and time periodic, whether or not the hydrodynamic interactions are included in 
the model. For short times, in case of the constant force, we show that the average drift is linear in 
time for long range hydrodynamic interactions and behaves ballistically or exponentially for local 
hydrodynamic interactions. Moreover, it can be opposite to the direction of external disturbance for 
some values of the model's parameters. When the force is time periodic, the effects are macroscopic: 
the system splits into two distinct spatial regions whose size is proportional to the value of the 
applied frequency. These two regions are characterized by different amplitudes and phase shifts in 
the response dynamics. 



I. INTRODUCTION 



Semiflexible and flexible polymeric chains [l|-_ 
branes [1, l8Hll| , moving interfaces |12h16J , growing sur- 
faces [lTj , single file systems [18[ are just few among the 
linear elastic systems whose time evolution is ruled by the 
stochastic equation 



8 



h(x,t) = / d d x'A(x-x") 



d z 



-h(x',t) + r)(x,t), 



dt v ' ' J v ' d\x' 

(1) 
termed as generalized elastic model [19|, |2f|- In 
its general formulation, the GEM (JXJ) is given for 
a D-dimensional stochastic field h defined in the d- 
dimcnsional infinite space x. The Gaussian random 
noise source satisfies the fluctuation-dissipation (FD) re- 
lation (r)j (x, t) r) k (x', t')) = 2k B TA {x - x') 8 3 k S(t - t') 
(j,k £ [1,D]), where A (r) = 1/ \r\ a corresponds to the 
hydrodynamic friction kernel whose the Fourier trans- 

( 47r )rf/2 r((d-q)/2) i ^a-d A] -,a-d 



form is A (q) 



I>/2) 



\q\ a - d = A\q\ a - d , if 



^2^" < a < d. The fractional derivative appearing in the 
right hand side of Eq.dTJ) is commonly defined as frac- 
as ._ _ /_Y72W 2 



tional Laplacian ( 



al3p 



(-vT 



211), and is ex- 



pressed in terms of its Fourier transform J-^ \ w=\* \ = 

— \qf [22J. A specific choice of the numerical values 
of the parameters characterizing Eq.([T]), namely D,d,a 
and z .leads to each one of the systems aforementioned 
(see [20] for a detailed description). However, a main 
partition between the systems obeying to Eq.([T]) can be 
done according to whether or not the hydrodynamic in- 
teractions can be considered long range. Indeed, if the 
hydrodynamics is only local, i.e. A (r) = 5 (f) , it must be 
formally set a — d and A = const in its Fourier transform 
expression A(q). For instance, setting D = 1, z = 4, 
a = 1 and d = 2 in ([1} corresponds to the stochastic 
equation for the height of a free fluid membrane float- 
ing in a solvent [1 S-EI, while for D = 3, z = 2, 



A (r) = 5 (r) Eq.([T]) governs the monomer's dynamics 
in a Rouse polymer [2j. 



A. Fractional Langevin equation scheme 

Tracing the derivation firstly obtained for the single file 
model [18j . recently we have shown [19| how to derive 
the stochastic equation governing the motion of a tracer 
or probe particle placed at the system's position x, that 
is the fractional Langevin equation (FLE) 



K + D^,h(x,t) = C(x,t), 



(2) 



with /3 — 
a 



[Wjj 



^- d and K+ = Tr^-isin^)— 4 „ 
d). Here the fractional Caputo derivative [2_ 
is defined as D p c <j){t) = j^gy f*^ dt' (t - t')' 13 ±(j> (t'j 
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dt' {t - t'Y 

(0 <; j3 < 1) and its Fourier transform is given by 
J~u} \ D^ip(t) > = {—iuY 4>(oj) . The fractional Gaussian 

noise C, (x, t) satisfies the second fluctuation-dissipation 
(FD) relation 



{a($,t)( k (x,t')) = k B T 



K^ 



T(l-f3)\t-t'\ 



7W 6j > k - 



(3) 



The FLE representation ^ is a mere change of wording 
compared to that furnished by (Q]). However, rephrasing 
the system's dynamics may result to be very helpful in 
the understanding of the physical picture subtenting the 
tracer's motion. This will be amply demonstrated by the 
following analysis. 



B. Generalized elastic model with localized 
potential 

In this article we generalize the outlined framework by 
deriving the tracer's FLE when a localized potential is 



applied to the probe in x* (tagged probe). Our starting 
point is the following GEM 



^-h(x,t) 



,ii 



Jd d x'A(x~x')x 
^h(x>,t) + F{h(x",t),t}6(x> 



- V (#> t) , 
(4) 
with the local applied force F. For such systems the 
tagged probe's FLE has been derived in two particular 
cases: membranes [l(| {z = 4, a = 1, d = 2), where 
the applied harmonic potential was introduced to mimic 
the action of an optical/magnetic tweezer; and single file 



systems [18( (z — 2 



= d = r 

constant, 



, where three types 
time-oscillating and 



of forces were analyzed: 
hookean. 

In Sec|TT]we draw the FLE for the untagged tracer (x ^ 
x*), besides the tagged one. Our analysis reveals differ- 
ent surprising regimes attained by the untagged probe 
in x, validating the Kubo fluctuation relations (KFR) 
in presence of a constant and time-periodic force in x*. 
Quoting Kubo [25|], these relations state general rela- 
tionships "between the response of a given system to an 
external disturbance and the internal fluctuation of the 
system in the absence of the disturbance". In Sec lIIII 
we study the case of an applied constant force. Indeed, 
while the tagged probe response is sublinear in time, the 
untagged probe's is always different at short times: it 
is linear for systems characterized by long range hydro- 
dynamic interactions, and it is ballistic (with a sign de- 
pending on the model's parameters) or exponential in 
case of local hydrodynamics. Asymptotically, the tagged 
and untagged probes undergo the same dynamical be- 
havior. In Sec lIVI the force is time periodic, in this case 
the macroscopic effects are persistent in time. Indeed the 
system splits into two macroregions whose size is defined 
by the value of the applied frequency loq: the responses 
of these regions markedly differs in amplitude and phase, 
allowing to make experimentally testable predictions on 
the viscoleastic properties of the system. In Appendix [X] 
we furnish the d-dimensional expression for the Fourier 
transform and its inverse. In Appendix [B] we report the 
theorem for the Laplace method for solving asymptotic 
integrals. In Appendix [Cl the formula for asymptotic so- 
lution of Fourier integrals is furnished. 



h (q, w) 



A? u ,{F{h(x*,t),t}}e 



v (q, w) 



\q\ d - a (-iu + A\qy /2 ) -iu + Affl' 2 ' 

where we made use of the short notation 7 = 2(z + 
a — d) . After multiplying both sides of the equation by 
K + (— ioj)° and, inverting the Fourier transform in space 
by means of Eq. (|A2[) . we get 



K+(-iu) h(x,uj) = 
T„ {F {h(x*,t),t}} 9 (|x - x*| , w) + C (x, u) , 



(6) 



where the non-Markovian Gaussian noise is introduced 
as [II 



C(£,w) 



+ OO 



dx"rj{x",Lj)<P{\x" -x\,uj). (7) 



The functions 9 and $ appearing in ([6]) and in (J7J are ex- 
pressed as 6 (|a?| ,w) = AZ a _ d / 2 (|x| ,u>) and $ (|x| , oS) — 
Zd/2(\x\,u), where 



2a (\S\,ui) = 



(8) 



(27r) d 

Their expression in time yields 



»+A\W 



K+\x\ 1 - d/2 n /3 f+ooji-.-A 

{ l l y/2 D p c f d\q\\q\ x 



replacing A = a — d/2 and A = d/2 respectively, and 
where 9(t) is the Heaviside step function, and Jd/2-1 de- 
notes the Bessel function of order d/2 — 1. We point out 
that both functions O and $ coincide in the case of local 
hydrodynamics (a = d and A = const). The physical 
picture behind the above mathematical derivation gets 
clear after inverting in time the equations ([6]) and ([7]): 



K + D^h(x,t) = 

X! oo dt'F{h(^,f) ) t'}e(|f-^|,t-f) + C(f,t), 

(10) 
with 



II. FRACTIONAL LANGEVIN EQUATION 



We start 
tion in the 

the Fourier 

r+00 



from Eq.Q and 
Fourier space, 
transform in space 



give its solu- 

We first define 

and time as 



h(q 7 u J )=f^d d xf^dth(x,t)e 



^ — i(q-x-ujt) 



, and intro- 



duce the notation for the time-Fourier transform of the 
force: T u {F{h(x*, t),t}} = /+~ dtF {h(x*,t),t} e*"*. 
We proceed to the FLE derivation for the untagged 
tracer in x. The solution of Eq. Q reads 



+00 



C{x,t)= I dx" dt'r){x",t')<&(\x' -x\,t-t'). 

J— OO J— OO 

(11) 
Thus Q (\x — x*\ ,t — t') can be seen as the propagator 
carrying the external perturbation, exerted at the point 
x* at time t', to the point x at time t. Likewise, the 
function $ (| x — x' | , t — i!) represents the propagator of 
the Brownian random source rj (x, t) from the point x! to 
the point x in the time elapsed between t' and t. 
We now turn to the derivation of the FLE for the tagged 
tracer. In this case it will be sufficient take the limit 



* in (J6J) , which corresponds to set \x\ = in 

d / 2 {\x\ ,ui) expression ([SJ. To this end we recall 



x — > x 
the I a _ 
that the Bessel function expansion for small argument 

-d/2 [--l 



is Jd/2-i(r) 



(f) 1 a ' Z [23 ), from which follows 
= 1 after straightforward pas- 



e(0,uj)=Al a _ d/2 (0,uj) 

sages. Substituting in ^ and inverting in time gives the 

FLE expression for the probe particle placed at x* : 



«"+Dgh(f*,t) = F{h(f*,t),t} + C(f*,*). (12) 

In the upcoming sections we analyze the cases of a con- 
stant and a time-periodic external force applied in x* . 



III. CONSTANT FORCE 



Let F n represents the force along one direction only, 
say F a = F 3 : F{h(x, t),t} — F a 9(t). We are interested 
in the average drifts {h(x,t)) Fo and (h(x*,t)) Fo where 
we dropped the index j. By averaging both Eq. (|10p and 
(TT21 and plugging in the definition ofO(|af— x*\,t — t') 
one has 



_ A\x— x*\ 



y.S*! 1 - 1 / 2 



Fofodt'i 



(h(x,t)) Fo - J^ym 

(h(x ,t))F — K+r(i+p) 



■iP. 



(13) 

Comparing the tracers' responses in (fT3"]) it is clear that 
the Einstein relation only holds for the tagged probe, i.e. 



(h(x*,t)) Fo = 



(S 2 Ht)) 
2k B T 



Fn 



(14) 



As a matter of fact, we recall that the mean square dis- 
placement in the absence of_ external force is given by 
(S 2 h(t)) = 2 ] 



k B T 



-t 13 LL9||. On the other hand the 

281. 



J K+V(l+p) 

untagged tracer fulfils the more general KFR [2_ 
and its time behavior presents very interesting features. 
Indeed, the analysis of Eq. (fT3"|) leads us to the following 
conclusions: (a) the drift (h (x, t)) Fo attains two different 
behaviors for times larger and smaller than the correla- 
tion time r = \x — f*| 7 ' /A; (b) for t < t, the response 
is dissimilar whether the hydrodynamic interactions are 
considered to be long range or local. 

• t <r. Long range hydrodynamic interactions.- In 
this case the integral over q appearing in (1131) can 
be performed asymptotically, the solution to the 
main order is 
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FIG. 1. (Color online) Constant force Fq. (a) Schematic 
representation of the response (|13p . Black (upper) line: Ein- 
stein relation (|14|l for the tagged probe. Grey line (second 
line from the top): average drift for the untagged tracer in 
presence of long ranged hydrodynamic interactions, dotted 
line represents the linear expression (|15[) at small times. Blue 
lines (third and fifth line from the top): {h(x,t))F for lo- 
cal hydrodynamics, upper (the third) and bottom (the fifth) 
curves stand for the positive and negative responses (|16[l for 
t < T di ff and 4m < z < Am + 2 and 4m + 2 < z < 4m + 4, 
respectively (dotted lines); the middle orange curve (fourth 
from the top) represents the drift for z = 2m, which is expo- 
nentially small oc e~ t ' Tdi f f at small times (see also panel (b)). 
Inset: representation of the phase Z7r/2 in ()16[). for z corre- 
sponding to the upper plane the drift has the same direction 
as _Fo while it is opposite in the bottom; exponentially small 
drift arises at z — 2m (orange (horizontal) solid circles), (b) 
Average probe's drift for the Edward- Wilkinson chain (D = 
1, z — 2,a — d — 1, and A = l/£, where £ is the damping). 
Simulations are carried out submitting the tagged tracer to a 
force Fo — 0.33 and detecting the average drift for different 
probes (solid green (middle) and magenta (bottom) curves). 
Dashed lines represent the theoretical expressions (|13|l : 

(h(x,t)) Fo =F [y/i/&Qe- v '-(\x-x*\/2)erfc(yj\ with 

y = \x — x*\ y/£/(4t). Statistical averages were taken over 
2000 realizations. Other simulation's parameters are £ = 1.0 
and feT = 1.0. The expression of the diffusion time r, that 
entails the transition from short to long time behavior, is given 
by t — \x — x*\ c;, which is r = 400 for the green (middle) 
curve and r = 900 for the magenta (bottom) curve. 



is evaluated in Appendix [C] The untagged probe's 
drift expression in Eq. (|13[) is given by 



(h{x,t)) Fo 



r(a/2) 2 a - d 



A 



r((d-a)/2) 7T rf / 2 \x-x*\' 



-F t. (15) 



{h{x,t)) Fo ~ 
^zT(m^) sin (^) r^Zf, (16) 



Local hydrodynamic interactions.- We put a = d 
(and hence 7 = 2z) in Eq . (fT3"|) . The integral over q 



for z ^ 2m, with m G N. The results ([15)) and 
P6p suggest that the untagged probe moves in av- 



erage as a free Brownian and ballistic particle re- 
spectively, under the influence of an external effec- 
tive force whose amplitude is inversely proportional 
to the distance \$ — x*\. However, the ballistic pic- 
ture is reductive in the case of local hydrodynam- 
ics. As a matter of fact for 2 + Am < z < A + Am 
the response of the probe is opposite to the exter- 
nal disturbance Fq, while for Am < z < 2 + Am 
they have the same sign. For z — 2m the re- 
sponse is slower than any power so that we expect 
(h (x, t)) Fo ex F tP +1 /T e~ T l l (up to numerical pref- 
actor in the exponential function). These surpris- 
ing and counterintuitive results are summarized in 
Figure Ufa) . 

To take an example, let us discuss the situa- 
tion of growing surfaces. The cases z = 2, 
z = 3 and z = A (d = 1) refer to different 
types of atomic diffusion on a crystalline surface 
[17j . The response of the step h(x,t) (the line 
boundary at which the surface changes height) 
to Fq grows in time exponentially in the case 
of attachment-detachment diffusion (z — 2): in 
this instance the complete solution of the first 
of Eq. (fT3"|) is achieved, yielding (h(x,t)) Fo — 

\/i7(<)e-« 2 - (\x - x*\ /2) erfc(y)] (Figure 



F, 



Up))), where y = | x - x* | VW*) , A = l/£ 
and £ is the viscous coefficient; its behavior 
at small times is found to be (h(x,t)) Fo ~ 

For the ter- 
gives a negative 



3/2 



-2 e -£|x-**| 2 /(4t) 



W/t) 

race diffusion (z = 3) Eq.([T 
drift, i.e (h(x,t)) Fa 



6 F (Aty 



Finally if 

T |x— X*| J 

■ 4, the so called periphery diffusion, one re- 
covers the exponential growth at short times, i.e 
ex F t 7 ^ \x - 



|-4 e _|x-x*| 4 /(At)_ 



• t >r. The integral in |<f| appearing in the first of 
Eqs. (fT3l can be solved by using the Laplace method 



(see Appendix flBJ) 1 ): (h(x,t)) Fo 



tP. 



K+T(l+B) 

Thus the Einstein relation (fT4")l is then regained 
only when t > r, namely when the correlation 
lenght £ = {At) 2 / 1 exceeds the distance \x — x*\. 
The transient violation of Einstein relation con- 
trasts with the second FD Eq.©, which is always 
fulfilled: the larger the distance \x — x*\, the longer 
the transient. 

We can summarize the results obtained in this section 
in the following compact form. The average drift is cast 

as 



z — d 



(h(x,t)) Fo =F \x-x*\ z - a f 



(17) 



The scaling function / [it] exhibits two distinct be- 
haviours whether m < 1 or ti > 1. From fli~5j) and ([i~6|) . 
it turns out that when u«l 



a— f" r (f) .. , 

„.d/2 r( d T a ) ' 

/M-.-i^, 8 in(^)r(|)r(^)^ a) ( 18 ) 

Hi) 



8+1 — — 

oc m mt e « 



for i) long range, ii) local (z 7^ 2m) and ra) local (2 = 
2m) hydrodynamic interactions, respectively. When u 3> 
1 we have invariably 



/[«] 



EC^U* 



2 d - 1 7r rf / 2 (^ - d)r (|) 



(19) 



Furthermore, it can be shown [29( that the following 
relation holds rigorously for both tagged and untagged 
tracers 



{h{x,t)) Fo = ^~ji F o, (20) 

which, in turns, encompasses the Einstein relation (fj"4")) 
and its generalization, namely the KFR. 



IV. TIME PERIODIC FORCE 

We now consider the force F {h(x, t),t} = Fq cos(uiot) 
and from (1101) and (TT21) we have 



K+D^h(x,t) = F Q $ie [/+~ dt'e'^^Q (\x - x*\ , t - t') 

C(f,i), 

K+D^h(x*,t) = F 3?e [e"^ '] +C(f*,t). 

( 21 ) 
In this case we study the complex mobilities or admit- 
tances /i (x,luq) and /1 (a?*, wo) which are defined through 
the linear response relations [25| 



(v(x,t)) Fo = 3fc [/i(x,w )^oe- Mot ] 



22) 



(v(x*,t)) Fo = 5Re ^(f*, Wo )F e-^*] . 
Both tracers fulfill the generalized Green-Kubo relation 

(v(x,uj )v(x*,ui' )) = 4ft8(ui + u' )k B T$le{n (x,ui )] , 

(23) 
where the complex mobilities follow from Eq. (|2"Tj) 
through the definition (|2"2"|) . 



/- n A\x-x*\ 1 ~ d/2 r+00 , 1 - I -jQ-d/2 

M(g»tt>o)= (27r) d/2 — Jo d k1l9l 

J^-tdgllf-^]) — g^ , 



(24) 



and the unperturbed velocity correlation function on the 
left hand side of Eq. (f2"3"f is obtained from Eq. © . In par- 
ticular, we can write the velocity autocorrelation func- 
tion in the frequency domain as (v (x*,u) V (a?*,u/)) = 



fcaT sm vy ' u 1 ~^4:Trd(u + uj'), and verify the relation 
(P5|) for the tagged probe: this is the standard (canoni- 
cal) formulation of the first FD relation [25| ■ 
We now analyze the low and high frequency behaviors of 
M(*,wo) = |/*(2,wo)|e-W' w °>. 

• u) t < 1. Changing variable (y — (A/loq) 2 ^' 1 \q\) 
in 5fte [/i (x, uiq)} and Qm [/i (x, wo)] and using the 
Bessels function's expansion for small arguments, 



one has jj, (x,ioo) 



,1-/3 



/& 



-i(l-/9)4 



• Wot > 1. Long range hydrodynamic interactions. - 
Performing the integral in q for 5fte \p (x, Wo)] and 
377i [/i (a?, wo)] we obtain for the response's ampli- 
tude 



|/z(x,u> )| 



r(a/2) 



2 a 



-rf 



.4 



T((d-a)/2) n d / 2 \x - S*\ c 



(25) 




ln|n(x,co )| 




(b) 



n |x-x*| 



while the phase is negligible, i.e tan ip (x, luo) ~ ^j_. 

Local hydrodynamic interactions.- We employ the 
same technique as for expression (|16[) achieving, for 

z ^ 2m (m € N), 



|/i(£,cj )| 



-l + d/2 



J r(f)r(^)|sin(f) 



|z-ir*| d 



(w t)" 



(26) 



and tan if {x, uj ) ~ —i^qt/ cos (4^) ■ When z = 2m, 
|/i (x, cjo)| is exponentially small and ip (x,u>q) ~ 
(1 — P)j — co7is/j(aj T) 1 / 2 : for instance the ex- 
act solution for z — 2,d — 1 reads /j(x,ujq) = 



[\$-x*\ e -<-{ 



\\i 



*l 



Let us now analyze in detail the physical scenario emerg- 
ing from the former analysis. The graphical rendering 
of the following discussion is presented in Fig 121 For a 
given frequency loq the system is divided into two spatial 

regions, \x - $*\ < \±) (I) and \x~x*\> (^J 

(II), characterized by very distinct dynamical phases. In 
case of long range hydrodynamics, the response of the 
system's portion closer to the tagged probe (I) shows a 
dependence of the amplitude oc cj and a phase shift 
(1 — /?)7r/2 with respect to the applied oscillatory force 
(already noticed in single-file systems 18]). On the 
other hand, in the outer region (II) the response's am- 
plitude decays as \x — x*\~ a , but almost no phase delay 
is displayed if compared to the external force. When the 
hydrodynamic interactions are local, although region (I) 
exhibits the same behaviour as in long range interacting 
systems, in region (II) the amplitude of the response is 
smaller and decays faster, namely oc \x — x*\ to^ if 

_,^_^*| 1/z 

z ^ 2m and oc e~* '"° for z = 2m; the phase shift 
instead is p>(x, ujq) sa ±-| for z ^ 2m, and it grows like 

—u! Q \x — x*\ if z = 2m. 




(0 



|x-x*| 



FIG. 2. (Color online) Time-periodic force Fq cos(ujot). (a) 
3D rendering of a membrane described by fTJ, under the ef- 
fect of an applied time-periodic force in x* (black arrow). 
h (x, i) represents the height of the fluctuating membrane on 
a 2-dimensional substrate (x — (x,y)). Regions I and II cor- 
respond respectively to the inner and outer region in which 
the membrane separates when the force is applied. The color 
code, red for region I and green for region II, has been drawn 
for the reader's convenience: increasing the frequency uio en- 
tails the shrinkage of the red region (I), (b) Schematic repre- 
sentation of the untagged response amplitude \fi(x,uJo)\ as a 
function of the distance \x — x*\: since is a schematic draw- 
ing no scale is needed on the :r-axis. In region I the universal 
behavior ljq~ /K + holds for any kind of hydrodynamic inter- 
actions. In region II the decay of the response's amplitude is 
oc \x — x*\~ a for long range hydrodynamic systems (Eq. (|25|) . 
grey (upper) solid line), oc \x — x*\~ z ~ for local hydrody- 
namic systems with z ^ 2m with raeff (Eq.(f26]), blue (mid- 
dle) solid line), and exponentially fast for local hydrodynamic 
systems with z — 2m (orange (bottom) line), (c) Schematic 
representation of the untagged response phase ip (x,uo) as a 
function of the distance \x — x*\. No scale is needed on the 
a;-axis. Region I: the system displays an universal phase de- 
lay for long range and local hydrodynamic interactions, i.e. 
tp (x,uio) = (1 — /3)f- Region II: For long range hydrodynamic 
systems the phase is absent, i.e. <p(x,ujo) ~ (grey (mid- 
dle) solid line); for local hydrodynamics the phase is approx- 
imately 7r/2 if 1 + 4m < z < 3 + 4777 with m 6 N (upper blue 
line), while it is approximately —tt/2 if 4m < z < 1 + 4m or 
3 + 4777 < z < 4 + 4m (bottom blue line); if the hydrodynamic 
interactions are local and z = 2m the phase shows a linear de- 
pendence on the distance \x — x*\, i.e. ip (x, cjo) ~ — \x — x*\ 
(orange (bright linear) solid line). 



CONCLUSIONS 



In this paper we derived the FLE representation of the 
tagged (x*) and untagged (x) probes' dynamics when 
a localized potential acts on x* . We demonstrated the 
validity of the KFR and generalized Green-Kubo rela- 
tion for both tracers, and the ensuing non-trivial physical 
regimes. This findings have important experimental and 
theoretical consequences. 

From the experimental point of view, the response to a 
constant force exerted on a position x* (implemented by 
an atomic force microscope by instance) can be detected 
experimentally within the domain of single-particle track- 
ing [30]. Indeed, the motion of the untagged tracer (a?), 
be an optical label, such as a gold or polystyrene bead, or 
a fluorescent tag, may provide a direct probe of the vis- 
coelastic properties of the system under study, as well as 
of its underlining elastic energy. The single-particle pos- 
sible responses are well schematized in Figure [TJ where 
the different time behaviors undergone by the average 
drift are displayed. In particular, we notice the surpris- 
ing effect for which the untagged tracer moves opposite 
to the external force for short times. 
On the other hand, our analysis provides a quantitive 
clear-cut description of the macroscopic observable ef- 
fects that a localized perturbation produces on an elastic 
system modeled by ([T]). As a matter of fact, a local oscil- 
lating field separates the systems in two regions, whose 
size can be tuned by tuning the amplitude of the charac- 
teristic frequency ujq. Moreover, we demonstrated that 
the behavior attained by these macroscopic domains is 
characterized by very different amplitudes and phases, 
according to the type of interaction and the values of 
the parameters which set the model ([T]). Figure [5] shows 
the measurable and testable predictions that our analysis 
enucleated. By instance, the readout of the effects of a 
local perturbation exerted by AFM, could be done using 
differential confocal microscopy to image the membrane 
ripples [3l[ . 

We believe that these findings have important direct ap- 
plications in the biosensors design and single-molecule 
manipulations. 

From the theoretical point of view we have shown that 
the FLEs (fH))) and (TT2"]) constitutes a powerful and com- 
prehensive dynamical representation of the motion of 
both tagged and untagged probes. Indeed, this stochas- 
tic equatio n, ly ing in the class of generalized Langevin 
equations [25], can be seen as a stochastic representa- 
tion of the general Kubo fluctuation relations with the 
single-probe random force £(a?, t) satisfying the second 
FD relation. 



Appendix A: Fourier transform of a d-dimensional 
isotropic function 

The d-dimensional Fourier transform of <f)(f) which is 
function only of its modulus |r|, i.e. <f>(f) = </>(|r|), is [321 ] 



(27r)d/2 \tf-W f o +°° d | r1 | r f / 2 Ja/^flsl |r1M|r1). 

(Al) 
Conversely, its inverse Fourier transform 0(|(?|) is given 
by 



c (i^ r >(ki) 



\r\ 1 - d ' 2 f +o 
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Appendix B: Laplace method for asymptotic 
integrals 

We hereby report theorem for the asymptotic solution 
of exponential integrals through Laplace method [33| . 
Consider the integral 



,6 

I(i) = / dqx(q) e 

J a 



p{q)t 



(Bl) 



If the following hypothesis are fulfilled 

i) p(q) > p(a) for any q £ (a, b) and the minimum of 
p(q) is appraoched only at a; 

ii) ^ and x(q) continuous functions in a neighbor- 
rod of a, except, possibly, at a; 

in) as q — » a + , (p(q) — p(a)) ~ P (q ~ a)^ and x(q) ~ 
Q (q — a) A_1 , where P, /i and A are positive constant and 
Q e R or e C; 

iv) I(t) is absolutely convergent throughout its range 
for all sufficiently large t; 
then the integral I(t) is 



I(t) 



0, 



-p(a)t 



M \fij (Pt) X/f * 



(B2) 



in the limit t — > oo. 

The integral in (1131) satisfies the previous hypothesis with 

a = 0, b = oo, P = A, Q = r/r1/r ,L~ < ,n- d / 2 , V = 7/2 



and X = a. 



r(d/2)\x-x*\ 1 



Appendix C: Asymptotic solution of Fourier 
integrals 



We hereby show how to derive Eq. ([16[) from the cor- 
responding general formulation of the untagged drift in 
Eq. tfTB")) . We first apply the change of variable y — 
(At) 2 / 1 \q\, achieving 



iv. 1^ x\\ A\x-x*\ 1 - d/2 x-, c-t ,,, 

(h(x,t)) Fo = (27r) l /2 — F J Q dt'x 
Of)^ I +O °dyy^J d/ ^ 1 (Xy)e-y\ / 



(CI) 



where A = \x — x*\ /{At') 2 / 1 is a large parameter. The 
integral over y can be evaluated by expanding the expo- 
nential for small arguments, i.e. e~ v ~ 1 — y z : the first 
term gives zero contribution, while the second is |26| 



ti 00 dvv*+ d '*J d/2 _ 1 (\y) 
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Alternatively, we can use the value of the improper inte- 
gral 



+00 



dyy v e 



u p -i\ v _ r(t/+l) iK (v+1) 



X v+i 



(C3) 



with v > —1, which can be obtained by the method of 



summation of improper integrals [34JJ. Recall that the 
Bessel function for d = 1 is J_i/2 (x) = W^-cosx, for 

rf = 3 is J1/2 (x) = w^sinx and for d — 2 Jo (x) ~ 

v /X C os(a;-f) for large as Ht 

The real and imaginary part of the mobility /^ (x, Wo) 



(C2) are obtained in the same way. 



ACKNOWLEDGMENTS 

A.C. and J.K. acknowledge the support of Marie Curie 
IIF programme, grant "LeFrac" . A.T. is deeply indebted 
with Tod Bertuzzi (TOY) , Gianluca Schneider-Faberi aka 
"ercapoccetta" and Andreas Gebardth for their great 
work on Figure [2j A.T. aknowlegdes Dario Villamaina, 
Andrea Puglisi, Angelo Vulpiani and the TNT group for 
illuminating discussions. 



[1] Doi M and Edwards SF, The Theory of Polymer Dynam- [20 

ics (Clarendon, Oxford, 1986). 
[2] Rouse PE, J. Chem. Phys. 21 1272 (1953). [21 

[3] Zimm BH, J. Chem. Phys. 24 269 (1956). 
[4] Granek R, J. Phys. II France 7 1761 (1997). 
[5] Farge E, Maggs AC,Macromol. 26 5041 (1993). [22 

[6] Caspi A, Elbaum M, Granek R, Lachish A, Zbaida D, [23 
Phys. Rev. Lett. 80 1106 (1998). [24 

[7] Amblard F, Maggs AC, Yurke B, Pargellis AN, Leibler 

S,Phys. Rev. Lett. 77 4470 (1996). [25 

[8] Freyssingeas E, Roux D, Nallet F, J. Phys. II France 7 [26 

913 (1997). 
[9] Heifer E, Harlepp S, Bourdieu L, Robert J, MacKintosh [27 
FC, Chatenay D, Phys. Rev. Lett. 85 457 (2000). 
[10] Granek R, Klafter J,Europhys. Lett. 56 15 (2001). 
[11] Zilman AG, Granek R, Chem. Phys. 284 195 (2002). [28 

[12] Edwards SF, Wilkinson DR, Proc. R. Soc. London A 381 

17 (1982) . [29 

[13] Joanny JF, de Gennes PG, J. Chem. Phys. 81 457 (1984). [30 
[14] Searson Rong Li PC, Sieradzki K, Phys. Rev. Lett. 74 

1395 (1995). 
[15] Krug J Scale Invanance, Interfaces and N on- Equilibrium [31 

Dynamics (Plenum, New York, 1995). 
[16] Krug J, Adv. Phys. 46 139 (1997). ' [32 

[17] Toroczkai Z, Williams ED, Phys. Today 52, No. 12 24 

(1998). [33 

[18] Lizana L, Ambjornsson T, Taloni A, Barkai E, Lomholt 

M, Phys. Rev. £81 051118 (2010). [34 

[19] Taloni A, Chechkin A, Klafter J, Phys. Rev. Lett. 104 
160602 (2010). 



Taloni A, Chechkin A, Klafter J, Phys. Rev. E 82 061104 
(2010). 

Samko SG, Kilbas AA, Marichev OI, Fractional Integrals 
and Derivatives, Theory and Applications (Gordon and 
Breach, Amsterdam, 1993). 
Saichev A, Zaslavsky G, Chaos 7 753 (1997). 
Caputo M, Geophys. J. R. Astr. Soc. 13 529 (1967). 
Podlubny I Fractional Differential Equations (Academic 
Press, New York, 1999). 
Kubo R, Rep. Progr. Phys. 29 255 (1966). 
Abramowitz M, Stegun I, Handbook of Mathematical 
Functions (Dover, New York, 1964). 

Marconi UMB, Puglisi A, Rondoni L, Vulpiani A 
Fluctuation-dissipation: response theory in statistical 
physics. Phys. Rep. 461 111 (2008). 
Villamaina D, Puglisi A, Vulpiani A, J. Stat. Mech. 
L10001 1 (2008). 

Taloni A, Chechkin A, Klafter J, to be published. 
Saxton MJ, Fundamental concepts in biophysics: Volume 
1. Handbook of modern biophysics, ed. T. Jue (Humana 
Press, New York, 2009). 

Lee CH, Tsai FC, Wang CC, Lee CH, Phys. Rev. Lett. 
103, 238101 (2009). 

Champeney DC, Fourier Transforms and Physical Appli- 
cations (Academic Press, London, 1937). 
Olver FWJ, Asymptotics and Special Functions (Aca- 
demic Press, New York, 1974). 

Hardy GH, Divergent Series (Clarendon Press, Oxford, 
1949). 



